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ABSTRACT 

We imaged five objects near the star forming clouds of Ophiuchus with the 
Keck Laser Guide Star AO system. We resolved Allers et al. (2006)'s #11 
(Oph 16222-2405) and #16 (Oph 16233-2402) into binary systems. The #11 
object is resolved into a 243 AU binary, the widest known for a very low mass 
(VLM) binary. The binary nature of #11 was discovered first by Allers (2005) 
and independently here during which we obtained the first spatially resolved 
R rsj 2000 near-infrared (J & K) spectra, mid-IR photometry, and orbital motion 
estimates. We estimate for 11A and 11B gravities (log(g)> 3.75), ages (5 ± 2 
Myr), luminosities (log(L/L Q ) =— 2.77±0.10 and — 2.96±0.10), and temperatures 
(T eff = 2375 ±175 and 2175 ±175 K). We find self-consistent DUSTY evolution- 
ary model (Chabrier et al. 2000) masses of l7~^Mj up and 14+ 6 Mj Mp , for 11A and 
11B respectively. Our masses are higher than those previously reported (13-15 
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Mj up and 7-8 Mj up ) by Jayawardhana & Ivanov (2006b). Hence, we find the 
system is unlikely a "planetary mass binary", (in agreement with Luhman et al. 
2007) but it has the second lowest mass and lowest binding energy of any known 
binary. Oph #11 and Oph #16 belong to a newly recognized population of wide 
(> 100 AU), young (< 10 Myr), roughly equal mass, VLM stellar and brown 
dwarf binaries. We deduce that ~ 6 ± 3% of young (< 10 Myr) VLM objects 
are in such wide systems. However, only 0.3 ± 0.1% of old field VLM objects are 
found in such wide systems. Thus, young, wide, VLM binary populations may 
be evaporating, due to stellar encounters in their natal clusters, leading to a field 
population depleted in wide VLM systems. 

Subject headings: instrumentation: adaptive optics — binaries: general — stars: 
evolution — stars: formation — stars: individual (2MASS J16222521-2405139 
and 2MASS J16233609-2402209) — Brown Dwarfs — extrasolar planets 

1. Introduction 

The system 2MASS1207334-393254 (sep~41 AU, 2M1207A~21M Jup , 2M1207b~5M Jup ) 
demonstrated that very low-mass objects could remain bound for at least 8 Myr (Chauvin 
et al. 2004; 2005c; Mamajek 2005; Song et al. 2006). A somewhat similar very wide (~241 
AU), very young, nearly equal mass, brown dwarf binary (2MASS Jl 101192-773238) has 
been discovered in Chameleon (Luhman 2004). While Billeres et al. (2005) discovered a 
wide (~ 200 AU), much older, field binary Denis J0551-44. They noted that "wide ultracool 
binaries are undoubtedly rare in the field". Indeed, of the ~69 field VLM binaries now 
known only Denis J0551-44 and Denis 2200-30 (~ 38 AU) have separations wider than 
~ 30AU (Close et al. 2003; Bouy et al. 2003; see also Burgasser et al. 2006; and Siegler et 
al. 2005 and references therein). 

Calculations of opacity-limited fragmentation in a turbulent 3-D medium yield minimum 
masses ~7 Mjup (e.g., Low & Lynden-Bell 1976; Boyd & Whitworth 2005, Bate 2005 and 
references therein). Thus, a binary system with individual masses in the vicinity of 7 times 
that of Jupiter is theoretically plausible; however, it should also be very rare since even a 
> 5% fraction of higher mass binary brown dwarfs is difficult to produce from star formation 
simulations (Bate et al. 2002; 2003). Jayawardhana & Ivanov (2006b) recently claimed 
discovery of just such a wide "planetary mass binary". We independently discovered this 
very same pair "Oph 11" in the course of our high-resolution survey of the Allers et al. 
(2006) Ophiuchus dark cloud low-mass objects (we note in passing that Allers (2005) was 
the first to resolve this system). 
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Allers et al. (2006) present a near-infrared (1.2, 1.6, and 2.2 /xm; hereafter NIR) imaging 
survey of several young star formation associations. By cross-correlation with the "Cores 
to Disks (c2d)" Spitzer [3.6], [4.5], [5.8], [8.0], & [24] pim legacy mid-IR survey (Evans et 
al. 2003) of these same star formation regions, they were able to identify 19 candidate very 
young, low-mass, objects each of which appears to have a significant mid-IR excess in the 
Spitzer dataset. 

In a high spatial resolution imaging survey of the Allers et al. (2006) Ophiuchus cloud 
sources, we found their sources #11 & #16 to be split into 1.7 — 1.9" binaries. In sections 
4.1-4.3 we prove that #11AB is a physical binary. We consider the binary nature of #16AB 
in Section 4.5. 

While Allers et al. (2006) estimate a mass of ~ 9 Jupiters for Oph 11B, based on their 
optical spectra and NIR photometry, Jayawardhana & Ivanov (2006b) estimate masses of 
13-15 Jupiters and 7-8 Jupiters for 11A and 11B, respectively, with the DUSTY models of 
Chabrier et al. (2000). We fully summarize all the recent results on Oph #11 in section 4.4. 

We will show, after an analysis of each component's luminosity, surface gravity, age, 
and T e ff, that the masses of 11A and 11B are likely just above the deuterium-burning limit 
(~ 13Mj up ). While we find it unlikely that Oph #11 is a planetary-mass binary, we do 
confirm that it is unusually wide, and of extremely low-mass, with likely the lowest binding 
energy of any known brown dwarf binary. Oph #16 has the 4 th lowest binding energy of any 
known VLM system. 

In section 4.6 we note that #11 and #16 join four other recently discovered young, 
wide, VLM binaries which, we argue, define a new class of "young, wide, VLM binaries". 
We find such wide VLM systems to be very rare in the field but ~15 times more common 
in very young clusters and associations. In section 4.6.3 we suggest that stellar encounters 
in these clusters have the potential to break- up/evaporate these binaries. 

2. OBSERVATIONS & REDUCTIONS 

2.1. Imaging Oph #11 and Oph #16 at J, H, Ks, & Ls 

We independently discovered the binary nature of the Oph 11 system with the Gemini 
NIRI camera (Hodapp et al. 2003) on July 1, 2006 (UT) in the J, H, & Ks bands (all other 
objects in our survey were observed with the Keck II LGS AO system). We then obtained 
Ls (~ 3.0/im) images with the NIRC camera on the Keck I telescope on July 7, 2006 (UT). 
In all cases seeing was excellent, yielding images with FWHM~0.3" in the Ks band. We 
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reduced these data in the standard manner (Close et al. 2002, 2003). All images were fully 
flat-fielded, sky and dark subtracted, bad pixel masked, aligned, and medianed. 

Figures 1 and 2 and Tables 1-3 present the measured and derived characteristics of the 
Oph 11 and 16 binaries. Since the northern component is brighter (in the NIR) we name it 
Oph 11A and the apparently cooler southern component Oph 11B. 

The system is also known as 2MASS J 16222522-2405138 and Jayawardhana & Ivanov 
(2006b) refer to it as Oph 162225-240515 (and Oph 1622). We will refer to this system 
as Oph #11 (and as Oph 1622-2405) since it was first mentioned in the widely available 
literature as #11 by Allers et al. (2006) to be an interesting low-mass object (in fact it 
was 11B that Allers was referring to; Allers et al. 2007). The designation of Oph 1622 is 
inappropriate since Allers et al.'s #10, #11, #12, and #13 objects would all also be "Oph 
1622" due to their similar RA. 

The night of July 1, 2006 (UT) was photometric and zeropoints from photometric stan- 
dards were used for calibration at J, H, and Ks. There are small (2MASS — Gemini = 
—0.12 ± 0.09, 0.16 ± 0.10, 0.06 ± 0.10 mag) differences between our calibration and the J, 
H, and Ks fluxes of #11 in the 2MASS point source catalog. However, the 1.9" separation 
of the binary likely leads to some errors in the 2MASS 4" aperture photometry, and so we 
adopt our values (which take the separations of the binary into account) throughout the 
present paper. 

We also analyzed the public Spitzer [3.6], [4.5], [5.8], [8.0] /im IRAC images of #11. 
Utilizing a custom IRAC PSF fitting program we were able to detect (for the first time) 
both components of the binary in all IRAC passbands (see, for example, Figure 3). In all 
cases a double PSF model of a 1.94" binary with PA=176° was a significantly better fit to 
the IRAC images than a single IRAC PSF. 

2.2. NIRSPEC J and H band Spectra of 11A and 11B 

Spatially resolved J- and fT-band spectra of the #11 system components were obtained 
on August 9, 2006 with NIRSPEC in low resolution mode on the Keck II telescope (Mclean 
et al. 1998, 2000). We used a 0"57x42" slit and the image rotator to maintain the ~180° po- 
sition angle with both components on the slit. The 0"57-wide slit corresponds to 3 pixels 
on the detector and a nominal resolving power of R~2000. J-band observations utilized 
the NIRSPEC-3 order-sorting filter (1.143—1.375 /im), and -fT-band observations used the 
K filter (1.996—2.382 /mi). For each filter a set of four frames was obtained in an ABBA 
nod pattern, nodding ~20" along the slit between the A and B positions. Integration times 
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for each exposure were 360 seconds in the J-band and 180 seconds in the i^-band. An AOV 
calibrator star, HIP 80224, was observed in each band for a 60 second integration at each nod 
position and at a similar airmass to facilitate the removal of telluric absorption features from 
the target spectra. Flat, dark, and arc lamp exposures were obtained with the calibration 
unit internal to NIRSPEC. 

The NIRSPEC data were reduced using the REDSPEC codefj] a package of IDL proce- 
dures developed specifically for the reduction of NIRSPEC spectra by S. Kim, L. Prato, and 
I. S. McLean. The reduction of low resolution NIRSPEC data with REDSPEC is described 
in detail in McLean et al. (2003), and summarized here. REDSPEC uses arc lamp frames 
to spatially and spectrally rectify the raw data and to determine a wavelength solution from 
observed neon and argon emission lines of known wavelengths. Following rectification, pairs 
of target nods are subtracted to remove background and divided by the rectified, dark- 
subtracted flat field. The target and calibrator spectra are extracted by summing over 9-11 
rows of data centered on the peak of the trace. For the Oph 11 observations, the peaks of 
the traces are separated by 18 pixels, resulting in minimum contamination between compo- 
nents. Intrinsic spectral features of the calibrator (Pa/5 at 1.282 /mi in the J-band and Br7 
at 2.166 /im in the A'-band) are removed by linear interpolation over the line. The calibrator 
spectrum is divided from the target spectrum, and the target spectrum is multiplied by a 
9770 K blackbody (the T e ff of HIP 80224) to restore the spectral slope. To produce the 
final spectrum, the extracted spectra from all four nods are averaged and normalized to the 
continuum level at a given wavelength. 

The actual resolution of the observations was estimated by fitting a Gaussian to several 
emission lines in the raw arc lamp frames, which resulted in R~1500 for J-band and R~1900 
for A-band. The signal-to-noise is calculated from the maximum peak-to-peak variation and 
estimated to be S/N > 100 for each spectrum. 

3. ANALYSIS 
3.1. Extinction 

The JHK colors from Table 1 are bluer in H-K and redder in J-H than the stellar locus 
of old M & L dwarfs (see Figure 4). Kirkpatrick et al (2006) discuss the peculiar shape of the 
H-band in the spectra of low gravity objects. These colors are not due to a highly (line of 
sight) extincted source since the visible spectrum of #11 was measured to be that of a young 
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M9 by Jayawardhana & Ivanov (2006a) and they (as do Allers et al. 2006; 2007) find A y =0 
is a good fit to the system. Moreover, our own comparison of the optical spectra of 11B from 
Jayawardhana & Ivanov (2006b) to a young L0 (2MASS J01415823-4633574; Kirkpatrick et 
al. 2006) also suggests the extinction towards 11 A and 11B must be very low. Consistent 
with a low Ay, we note there is little 24 pm IR cirrus in the immediate area around #11. 
This is not surprising as #11 is located (dist ~ 0.5° ~ lpc) near the edge of the 13 CO core 
of the p Oph cloud, although #11 (and #16) are still inside the p Oph survey area (and the 
13 CO cloud) of Wilking et al. (2005). 



3.2. J and K Band Spectra of 11B 

From the mainly gravity independent, yet temperature sensitive, K-band CO lines (see 
for example Gorlova et al. 2003) we can measure T e ff and spectral type for 11A and 11B. 
However, this first requires some calibrated, young, cool brown dwarfs for comparison. There 
does not yet exist a standard set of very young, late M/early L spectral standards. However, 
the young brown dwarf 2MASS J01415823-4633574 (hereafter 2M0141) has been studied with 
relatively high resolution spectra from 0.5-2.5 pm by Kirkpatrick et al. (2006). They find 
that 2M0141 is best fit by an L0 spectral type and by T eff = 2000±100K and log(g)=4.0±0.5 
from detailed optical to NIR spectral synthesis fits. Hence, this is an obvious brown dwarf to 
compare our spectra to (see Figures 5-8). In addition, McGovern et al. (2004) obtained some 
J-band spectra of a series of young brown dwarfs with a similar instrumental configuration 
(hence spectral resolution) of NIRSPEC as we did. Hence, we will also compare our J-band 
spectra to theirs. 



3.2.1. The Temperature of 11B 

In Figure 6 we compare 11B to 2M0141. The excellent match of HB's CO, Nal, Call, 
and pseudo-continuum to the K-band spectrum of 2M0141 suggests an T e ff ~ 2000K and 
~ 5Myr age (log(g)=4). 

When we compare 11B in the J-band to 2M0141 (middle trace of Figure 8), 11B appears 
slightly too hot for a good fit. In fact, a better fit to the J-band of 11B is the young (5 ± 2 
Myr) somewhat hotter (M9) brown dwarf o Ori 51 (McGovern et al. 2004). Hence, it 
appears that in the J-band 11B appears closer to that of an M9. Uncertainty at the level of 
one subclass for young brown dwarfs is not unusual. Since the T e ff for an M9 is 2300K from 
the temperature scale of Martin et al. (1999), 2400K from that of Luhman et al. (1999), 
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and 2400K from Golomoski et al. (2004) we adopt a value of 2350K for M9. So it appears 
11B could be as hot as T e ff = 2350K from its J-band spectrum (but as cool as 2000K in 
the K-band). Therefore, we adopt T e ff = 2175 ± 175 (M9.5 ± 0.5) as a reasonable match 
to the J & K band spectra of 11B compared to other young brown dwarfs. We note this is 
consistent with the M9-L0 spectral type found by Jayawardhana & Ivanov (2006b) from the 
visible spectrum of 11B. 

3.2.2. The Surface Gravity and Age of 11B 

In both the J and K spectra, log(g) appears consistent with ~ 4.0 and age 5 ± 2Myr. 
Younger ages seem to be precluded for 11B due to the poor fit of 11B to the very young 
(1-2 Myr) KPNO Tau-4 brown dwarf. Since the M9.5 spectral type of KPNO Tau-4 is 
consistent with the M9-L0 of 11B, we reason that the very poor fit of 11B is primarily due 
to significantly higher gravity in 11B compared to KPNO Tau-4. Indeed it is well known 
that the J-band has a host of gravity sensitive features (like the KI doublets and its entire 
pseudo-continuum; McGovern et al. 2004), hence we can reconcile the poor fit of 11B to 
KPNO Tau-4 and the good fit to a Ori 51 by adopting an age of 5 ±2 (hence surface gravities 
of log(g)> 3.75 according to the DUSTY models) for the Oph 11 system. These ages are 
significantly higher than the 1 Myr assumed by Jayawardhana & Ivanov (2006b). 

However, older ages of ~ 30Myr are very unlikely since the strength of the Nal and Cal 
lines in 11B (Fig. 5 & 6) are much weaker than that of GSC 8047 (a known 30 Myr old M9.5 
Tuc association member; Chauvin et al. 2005b). 

3.2.3. The Temperature and Age of 11 A 

We find 11A is consistently ~200K hotter than 11B from detailed spectral synthesis 
modeling of the J-band and K-band spectra (section 3.3.1). Moreover, M8.5-M9.5 spectra of 
age 5 ± 2 Myr should be a reasonable fit to llA's slightly hotter spectra than the standards 
of Figs 5 & 7. Hence, a T eff ~ 2375 ± 175K (or M9 ± 0.5) is adopted for 11 A. Again our 
NIR spectral type is consistent with the M9 spectral type found by Jayawardhana & Ivanov 
(2006b) from llA's visible spectrum. Although it is somewhat hard to disentangle the effects 
of temperature and surface gravity, the similarity of llA's spectra to that of 11B suggests 
that the age of 11 A is consistent with 5 ± 2 Myr. 
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3.3. Synthetic Spectral Fits 

3.3.1. Temperatures 

We consider how our adopted ~5 Myr age and M9 and M9.5 spectral type estimates 
compare to synthetic spectra with log(g)~3.95 (consistent with the DUSTY model's 5 Myr 
isochrone) and T e //=2375K for 11A and T e //=2175K for 11B. In Figure 9 we compare our 
observed K-band spectra to synthetic spectra computed using up-to-date PHOENIX DUSTY 
atmosphere models (Hauschildt et al. in prep). The synthetic spectra do a reasonable job 
fitting Nal and Cal, whereas the CO is a little stronger in the synthetic spectra, suggesting 
higher temperatures for 11A and 11B. However, we note that at these gravities some over 
estimation of the temperatures from the synthetic CO lines were also seen in the similar case 
of 2M0141 (Kirkpatrick et al. 2006). 

If we let T e ff and log(g) be free parameters, the very best fit to our observed spectra 
with the synthetic models suggests that AT e // ~ 200K between 11A and 11B. Hence, even 
though there is some uncertainty in the absolute temperatures derived from these fits, we 
can have some confidence in a AT e jj = 200K between 11A and 11B. 

3.3.2. Surface Gravities 

In figure 10 the J-band is compared to synthetic spectra. In the more gravity sensitive 
J-band we find our adopted log(g)=3.95 for Oph 11 is slightly too low for a good fit to 
the observed pseudo-continuum. A slightly better fit is obtained with log(g)=4.25 ± 0.50. 
However, the adopted T e jf = 2375K and T e ff = 2175K for 11A and 11B combined with 
log(g)=3.95 produce good agreement on the strength of the KI doublets. Hence, it appears 
that our adopted age of ~5 Myr (hence log(g)~3.95) for 11A and 11B is reasonable. 

In summary, the optimal synthetic fits suggest somewhat higher gravities (log(g)=4.25± 
0.50) but confirm our AT e jj = 200K between 11A and 11B. Our synthetic spectra preclude 
log(g)< 3.75 and so ages < 2 Myr are very unlikely for Oph 11. 

3.3.3. The SED of 11 A and 11B and Their Thermal IR Excess 

Longward of 3 fim the SED of 11A and 11B (Table 1; Fig 11) indicates substantial 
excess emission compared to our synthetic M9 and M9.5 (log(g)=3.95) spectra for 11A and 
11B. This suggests that both 11A and 11B have circumstellar dust disks. In particular, HB's 
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excess is very strong. The discovery by Jayawardhana & Ivanov (2006b) that both 11A and 
11B have strong (and broad) Ha emission suggests active accretion may be still occurring 
around both of these objects. While some 50 ±12% of young brown dwarfs have an IR excess 
suggesting circumstellar disks, only ~ 16 ± 6% are estimated to be actively accreting (Bouy 
et al. 2006b and references within). Oph #11 is the first brown dwarf binary where there is 
strong evidence that both components are actively accreting. Similar conclusions were drawn 
for the young, single, very low mass, brown dwarf Cha 110913-773444, which also has an 
IR-excess (Luhman et al. 2005). As noted in Section 3.1, extinction toward Oph 11 appears 
to be small; thus, the two dusty disks are not close to edge on. 



3.4. Luminosity and Mass of 11A and 11B 

With an Ophiuchus distance of 125 ± 25 pc (de Geus et al. 1989), an observed Ks = 
13.92±0.05 mag, a K-band bolometric correction of — 3.20±0.15 mag (Golomoski et al. 2004; 
appropriate for M9±0.5), and noting that for an M9 star Ks — K ~ 0.03 mag (Daemgen et 
al. 2007), we derive a photospheric luminosity for 11A of log(L/L Q ) = —2.77 ±0.10. This 
luminosity excludes any contribution from the excess IR emission. Similarly, we find with 
BC K = -3.10 ± 0.15 mag for an M9.5 ± 0.5 that log(L/L Q ) = -2.96 ± 0.10 for 11B. 

Notwithstanding some uncertainty in theoretical evolutionary tracks at such young ages, 
and low masses, it is reassuring that the errors between the DUSTY tracks of Chabrier et al. 
(2000) and the handful of dynamical mass calibrated systems are not large when accurate 
spectra and K band photometry are known (Luhman & Potter (2006); Close et al 2007; 
Stassun et al. 2006). Hence we use our system mass and the DUSTY isochrones to estimate 
masses for 11A and 11B. From the HR diagram in Figure 12 we see that both 11A and 11B 
fall close to the 5 Myr isochrone predicted by the DUSTY models. Over the suggested age 
range of 3 — 7 Myr we find estimated masses from Figure 13 of 13 — 21M Jup and 10 — 2QM Jup 
for 11 A and 11B, respectively. We are not the first to caution that these masses, based on 
theoretical isochrones at very young ages, may have unknown systematic errors. As well 
the range of NIR photometry from Allers et al. (2006), Jayawardhana & Ivanov (2006b), 
2MASS, and our work suggests there may be some variability in the NIR photometry (which 
would not be surprising for a system so young). 



3.4-1- Our Masses Compared to those of Jayawardhana & Ivanov 



As is clear from Table 1, our A J, AH, and AKs values are close to that found by 
Jayawardhana & Ivanov (2006b). Hence, we should derive similar masses based on the 
above technique which was similar to theirs. However, our absolute photometric calibration 
is closer to that of the 2MASS point source catalog. In particular, our integrated Ks flux 
is just 6% brighter than the 2MASS value, whereas Jayawardhana & Ivanov (2006b) find 
values ~30% brighter. Their significantly brighter K-band gave them closer agreement to 
the more luminous 1 Myr DUSTY isochrone, while we are closer to the 5 Myr isochrone (see 
Fig 13). Whereas, Allers (2005) and Allers et al. (2006) find Ks=14.03 ± 03 mag for 11B, 
some 40% brighter than we do in this study. Hence, deriving masses primarily based on NIR 
luminosity for #11 may be problematic for two reasons: first, there may be an significant 
underestimation of mass when one uses J and H fluxes and the dusty models (Close et al. 
2005; 2006); secondly, #11 may be a young variable - so even the use of the more reliable 
K-band luminosity (Close et al. 2005; 2007) may be misleading. Hence we will also estimate 
a luminosity (and distance) independent mass based on our log(g) and age estimates in the 
next section. 



3-4-2. Masses from log(g) andT e ff 

It is reasonable to assume that the age of 11 A and 11B must be < 8 Myr since there is 
a strong IR-excess for 11B, strong Ha emission (Jayawardhana & Ivanov (2006b)), and both 
11A and 11B are a reasonable fit to the a Ori 51 brown dwarf whose age < 8 Myr (McGovern 
et al. 2004). However, our fits to the J-band of young, late M brown dwarf standards suggest 
that the ages must also be > 2Myr (witness the very poor fit of KPNO-Tau4). Moreover, 
the synthetic spectral fits find log(g)> 3.75 also precluding ages < 2 Myr for either 11A or 
11B. From our age, surface gravity, and temperature ranges we find masses of l7ttM Jup and 
lA~^lMj U p are appropriate (see Fig. 13). Since these masses avoid some of the pitfalls of 
a luminosity based mass estimate we adopt these mean (~ 5 Myr) values, but we increase 
the uncertainty ranges to be consistent with our HR diagram masses of 13 — 21Mj up and 
10 — 20M Jup for 11A and 11B. Hence our final adopted masses are V7+\M Jup and \^%M Jup . 
Since we do not include (the currently unknown) systematic errors of these models, our mass 
uncertainties are likely underestimates. These masses are both slightly above the ~13 Mj up 
deuterium limit, and above the 14-15 Mj up and 7-8 Mj up masses found by Jayawardhana 
k Ivanov (2006b). 
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4. DISCUSSION 

4.1. Are 11 A and 11B Members of the Ophiuchus Cloud Complex? 

Ophiuchus 11B possesses a strong mid-IR (5 & 8 pm) excess, and 11 A a weaker one (see 
Figure 11). Mid-IR excesses are common only among stars younger than ~8 Myr (Mamajek 
et al. 2004; Silverstone et al. 2006; Rhee et al 2006 and references therein). The strong 
Ha emission measured by Jayawardhana & Ivanov (2006b) for 11A and 11B suggests #11 
is actively accreting and associated with the molecular cloud. Indeed, Wilking et al. (2005) 
find active accreters (CTTSs) ranging in age from ~ 0.3 — 10 Myr in the r ~ lpc area around 
the p Oph cloud (including the location of Ophll & Ophl6). Our derived 5 ± 2Myr age of 
#11 is consistent with Oph ll's ~ 0.5°(~ lpc) distance from the p Ophiuchus cloud core 
where active star formation is taking place. Hence, it is very likely that 11A and 11B are 
members of the Ophiuchus/Upper Soc sub group of the Sco-Cen OB association. 

Based on the very low extinction (Ay ~ mag) to Oph#ll (despite being in the line 
of sight to the 13 CO cloud), we suggest that Oph#ll is at a distance of ~ 125 ± 25pc. In 
other words, Oph #11 is likely slightly in the front of the ~150 pc "halo" of ~ 3 Myr objects 
around the p Oph cloud core, which typically have A v > 1.5 mag (Wilking et al 2005). To 
be conservative we adopt distances of 100-150 pc for the system. 

4.2. Do 11A and 11B Form a Common Proper Motion Pair? 

We consider whether 11A might be a foreground dwarf; statistically its observed proper 
motion is smaller than the expected proper motion of a foreground dwarf. We can use 
archival POSS II images to determine if the pair's orientation on the sky is changing (due 
to 11 A being a relatively fast moving foreground object). The very red optical colors of 11B 
make it impossible to detect in the POSS II "IR" (RG715) images, only a very faint object at 
the current position of 11A is observed, while neither is seen in the POSS II "blue" (GG395 
filter). However, both 11A and 11B appear in the RG610 POSS II "red" filter images (Fig, 
14) due to their strong Ha emission (Jayawardhana & Ivanov 2006b). In the 1993.33 epoch 
red POSS II images the binary separation is 2.05 ± 0.20" at PA = 174 ± 6°) which may be 
compared to our images (Figs. 1 & 2; separation^ 1.943 ± 0.022"; PA=176.2 ± 0.5°). 

If #11A were a field M9 dwarf with Ks=13.92 mag, then it should be at ~ 40 ± 5 pc. 
Adopting a velocity of V tan ~ 30 km/s (the approximate median of late M field dwarfs; Gizis 
et al. 2000), would then imply a typical displacement of ~ 2" between 11A and 11B in the 
epoch 1993.33 DSS images. Since this conflicts (at the ~ 10a level) with what is observed, 
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we conclude that 11A is not a foreground field dwarf. Hence, 11A and 11B very likely form 
a common proper motion pair. 

4.3. Are 11A and 11B Physically Bound? 

Comparison of the epoch 1993.33 and current images shows that the differential dis- 
placement of 11A and 11B over the last 13.22 years is only ~ 3 ± 5 km/s, consistent with 
expectations for a physical system this wide (estimated period (2 ± l)xl0 4 yr; Table 3). 
However, since all members of the Ophiuchus association will have similar proper motions, 
it is, perhaps, possible that 11A and 11B are well separated along our line of sight but, by 
chance, close together in the plane of the sky. 

We can estimate how likely it is to observe two very cool (late M/early L) objects within 
2" of each other in a cloud like Ophiuchus. The density of such objects can be estimated for 
the survey of Allers et al. (2006); in 1700 sq. arcmin they found two ~ M8-M9 objects (#11 
and #12). Follow-up optical spectra by Jayawardhana & Ivanov (2006a) showed that (in the 
optical) #12 appears to be a background quasar. We in fact find #12 to be a 1.3" "optical 
binary" consisting of a faint z=2 QSO "12B" and a slightly brighter G-giant "12A"- from 
Gemini GNIRS IFU observations). Hence, one would expect in Ophiuchus to find one late 
M/early L object in every 1700 sq. arcmin. Therefore the odds of a chance alignment of two 
such objects within 2" is just ~ 2 x 10~ 6 . Thus, most probably, 11A and 11B are physically 
bound. 

Are they still a physically bound pair today? If they are now drifting apart at, say, 
the ~0.5km/s escape velocity of the system, then the current 237 AU projected separation 
suggests that they must have been bound until just ~2xl0 3 years ago. In other words, to 
explain the close proximity of the pair they must have been bound for > 99% of their lives 
and have only just become unbound, a very unlikely scenario. It is, however, possible that 
the Oph 11 system will become unbound in the future (see Section 4.6.3.). 

4.4. Our Masses and Spectral Types Compared to those of Other Recent 

Studies 

Several groups have recently reported spectral types and masses for Oph 11AB, Allers 
et al. (2007) find spectral types of M7 ± 1 and M8 ± 1 for 11A and 11B (from the low 
resolution R~300 NIR spectra of Allers 2005). They derived ages of ~ 20 Myr and masses 
of 65 and 35 Mj up from dusty models. A more recent work by the same group (Luhman et 
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al. 2007) find types of M7.25 and M8.75 and masses ~ 59M Jup and ~ 21M Jup for 11A and 
11B and, like us, they derive an age of ~ 5 Myr for the system. 

These two studies derive earlier spectral types than our M9±0.5 and M9.5±0.5 for 11A 
and 11B. Part of the disagreement is due to a lack of "gold-standard" spectral templates at 
~ 5 Myr ages. Hence, both of these studies use/compose templates of somewhat younger 
(~ 1 — 3 Myr) ages (and with different extinction corrections) which can lead to additional 
uncertainty distinguishing between age and temperature effects (for example: cooler older 
objects have similar alkali line strengths as younger hotter objects). Moreover, there are 
systematic errors in the temperature scale, adding further uncertainty to the final T e ff. 

Very recently Brandeker et al. (2006b) have also estimated higher masses than Jayaward- 
hana & Ivanov (2006b) from new NIR spectra. They confirm our spectral types and tem- 
peratures (2350 ± 150 K and 2100 ± 100 K for 11A and 11B; compared to our 2375 and 2175 
K). They derive these temperatures by a pure comparison to the DUSTY models of Baraffe 
et al. (2002) and the latest brown dwarf models of Burrows et al. (2006). Over an age range 
of 1-10 Myr Brandeker et al. (2007) find masses of 13t 8 4 M Jup and l0± 5 4 M Jup for 11 A and 
11B (slightly lower, yet consistent, with our 17+\M Jup and \\^%M Jup masses). 

In summary, 11A is likely above the deuterium burning limit (~ 13Mj up ) and the system 
age is likely closer to 5 Myr than 1 Myr. However, there is still significant disagreement as 
to the exact spectral type of 11A with a range of M7.25-M9 (and dusty model masses of 
~ 58 — 13M Jup ). In the case of 11B there is better agreement (M8.75-M9.5; and masses 
of 21 — lOMjup). All recent studies agree that Oph 11AB is best described as a young (~ 5 
Myr), very wide (~ 240 AU), low mass (M tot ~ 23 — 80Mj up ) brown dwarf binary. Moreover, 
the evidence presented here indicates that the system is also a bona-fide bound binary with 
both members likely possessing their own accretion disks. 

4.5. The Oph #16 System 

We resolved Allers' object #16 into two components at JHKs (Fig.l & Table 2). In a 
similar manner as for Oph #11 we were able to estimate Spitzer [3.6], [4.5], [5.8], and [8.0] 
fluxes for both components (Table 2). 

4-5.1. Extinction, Colors, and Spectral Types for Oph 16 A and 16B 

Since Jayawardhana & Ivanov (2006a) did not obtain optical spectra of #16 it is some- 
what uncertain what the true spectral types for 16A and 16B are. We note, however, that 
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16A and 16B appear significantly more luminous and bluer (warmer) than 11A or 11B. 
Moreover, #16 is just 0.3° (~ 0.6pc) from the Oph cloud core (D= 125 ± 25 pc), and well 
inside the main dark cloud containing the L1688 core. Hence, Oph 16 would likely have an 
age of ~ 1 Myr. If this is true than we need an extinction of Ay ~ 3 ± 1 to have these 
objects fall near the 1 Myr isochrone of p Oph. This is a reasonable extinction since Allers 
et al. (2006) estimated Ay = 2 for this object. 

After correcting our NIR photometry for the Ay = 3 ± 1 extinction towards #16 we 
find (see Table 2) for 16 A J-Ks=0.99 ± 0.31 which is consistent with M2-M8 spectral types 
(Kirkpatrick et al. 1999). Hence, we adopt its spectral type as M5.5 ± 3. Future (spatially 
resolved) spectroscopic observations are required to obtain better estimates for the spectral 
types of 16A and 16B |§. However, we can estimate temperatures for 16A and 16B from our 
dereddened NIR colors: For 16A T eff ~ 3000 ± 300K and for 16B T eff ~ 2925 ± 300K 
(Golomoski et al. 2004). In the same manner as #11, we find (with BCk = —2.95) 
log(L/L Q ) = -1.18 ± 0.10 for 16 A and log(L/L Q ) = -1.47 ± 0.10 for 16B (where BC K = 
—3.00; Golomoski et al. 2004). The luminosity errors of #16 take into account our 1 mag 
of uncertainty in the extinction correction. 



4-5.2. Are 16A and 16B members of the Cloud Complex? 

From the location of 16A and 16B in the HR diagram of Figure 12 it is likely that both 
components are members of the Ophiuchus cloud (consistent with their ~0.6 pc distance from 
the core). From the ~ 1 Myr isochrone (in Figure 12) we estimate masses of ~ 100Mj up for 
16A (likely a VLM star) and ~ 72Mj up for 16B (likely a high mass brown dwarf or VLM 
star). Allers et al. (2006) estimated a mass for #16AB of ~ 110Mj Mp . Spatially resolved 
spectra in the NIR will be needed to better constrain the ages of 16A and 16B. However, the 
strong mid-IR fluxes of both 16A and 16B in Table 2 suggests that 16A and 16B are both 
young VLM members of the Ophiuchus cloud. 



4.5.3. Are 16A and 16B Bound? 

We believe that 16A and 16B form a common proper motion pair since the orientation 
of the binary today (separation=1.696± 0.005", PA=218.29±1.00; Fig. 1) is consistent with 
its orientation in 1982.66 (separation=1.87 ± 0.2"; PA=216 ± 10°; Fig. 14). This is exactly 
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what we would expect for a system with a very long period of ~ 8xl0 3 years (see Table 3). 

Is it plausible that 16AB could be a foreground VLM binary? We estimate (in the 
same manner as Close et al. 2003) such a hypothetical foreground system would have a 
photometric distance of ~ lOpc to match 16A's brightness and very red (undereddened) 
magnitudes (Ks=10.56 and J-Ks=1.41). Such a VLM binary would have a period of ~ 170 
years and so in the 23.89 years between Fig. 1 and Fig. 14 one would expect ~ 50° of rotation. 
However, over the last 23.89 years the system has only changed by a APA = 1.9 ± 10.0° 
which shows that Oph 16 is not a foreground VLM system. Moreover, the IR-excess of 
16AB cannot be explained by a foreground field (old) system. In the same manner as for 
Oph #11 we estimate that the odds of finding two unrelated, nearly equal magnitude, low- 
mass, Ophiuchus members within 2" of each other is less than ~ 10~ 5 . Hence, we conclude 
that Oph #16AB is a newly discovered, wide (~212 AU), low mass (M tot ai ~ O.17M ) VLM 
binary system in the Ophiuchus cloud complex. 

4.6. Ophiuchus #11 and #16 Compared with Other VLM Binaries 

Figures 15 & 16 show how Oph #11 and Oph #16 compare with other known binary 
systems. Oph #11 is the second least massive binary known (with M tot ~ O.O31M ; just 
slightly more massive than 2M1207) and has the lowest binding energy of any brown-dwarf 
binary. All currently known low binding energy systems are very young (< 10 Myr; plotted 
as open circles in Figures 15 & 16). Hence, one might predict that systems like Oph #11 
will tend to become unbound as they age. 

4-6.1. How Common Are Very Wide Young VLM Binaries? 

It is interesting to note that very wide, low-mass, binaries may not be that rare when 
very young. We observed four (not counting #12) of the Allers et al. (2006) Oph sample 
(their #8, #10, #16) and the two M8 objects (GY 264 & GY3) of Wilking et al. (2005) 
with Keck II LGS AO and found #11 (with NIRI) and #16 (~ 33 ± 23% of the sample) 
were binaries, and 100% of our binaries were wide. 

We caution that this is a very small sample, dominated by small number statistics. 
Also no correction for the Malmquist bias of the Allers et al. (2006) sample has been made. 
Nevertheless, a significant fraction of the Allers et al. (2006) Ophiuchus sample are young, 
wide, VLM binaries. Moreover, Bouy et al. (2006) find in their NGS AO survey of slightly 
more massive binaries in the nearby Upper Sco region that 3/9 of their binaries were wide 
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VLM systems. One of these binaries, Denis PJ161833. 2-251750. 4AB, is spectroscopically 
confirmed (Luhman 2005) and another, USCO 1600AB, is likely also bound (Bouy et al. 
2006). In contrast HST/ACS surveys of Upper Sco (5 Myr; Kraus et al. 2005) found no 
young wide VLM binaries. Bouy et al. estimate (including the non-detection of Kraus 
et al. 2005) that (3/12) ~ 25% of VLM binaries in Upper Sco are wide. In summary, 
Bouy et al. find two solid wide (> 100 AU) VLM binaries out of 40 VLM objects; hence 
fvLM widehvoung ~ 5±2% in the Upper Sco OB association. 

The (> 4 AU) spatial resolution survey of the lower density Taurus association of Kraus 
et al. (2006) with HST/ACS found 2/22 VLM objects were binary, but none were wider than 
~ 6 AU. Combining our (2/6) result and the (0/22) result by Kraus et al. (2006) suggests 
that 2/28 VLM objects are wide binaries. Hence, fvLM widekyoung ~ 7 ± 5% in Taurus and 
Ophiuchus. This is very consistent with the 5^% found in the Upper Sco cluster. Merging 
all the surveys indicates 4/68 young (< 10 Myr) VLM objects that were observed at high 
resolution (sep > 5 AU) were found to be wide (> 100 AU) binaries. Therefore, we adopt 
fvLM wideSzyoung ~ 6 ± 3%. While even larger surveys with LGS AO and HST will be required 
to confirm these frequencies, it is clear that a significant (~ 6%) fraction of young VLM 
objects are formed in wide (> 100 AU) binaries. 



4-6.2. Is There a Dearth of Wide VLM Binaries in the Field? 

In this section we try to estimate the frequency of wide VLM binaries in the field 
(fvLM wideSco i d ) based on published surveys. After correcting for the selection effect that field 
objects are fainter than young objects of a given mass, we can derive the depletion of old 
field wide systems to that of young systems. We define X evap =f VLMwidejkyoung /f VLMwidelkold , 
as the fraction of wide VLM binaries that have "evaporated" as these binaries age. 

Of the 69 field (old; 0.5-10 Gyr) VLM (M tot < O.2M ) binaries currently known (VLM 
Binary Archive § ) only Denis 0551 is wider than 100 AU. Hence only ~ 1.4% of known field 
VLM binaries are wide. All searches to find wide VLM binaries (Close et al. 2003; Gizis et 
al. 2003; Bouy et al. 2003; Burgasser et al. 2003; Siegler et al. 2005; Allen et al. 2007; Reid 
et al. 2006; Burgasser et al. 2006 and references within) have all failed to find a single wide 
(> 50 AU) VLM field binary despite being sensitive to sep > 2 AU and masses> 13Mj up . 
The 132 VLM (M7-L8) study of Allen et al. (2007) alone suggests the frequency of wide 
VLM systems must be fvLM widekold < 2.3% in the field at the 95% confidence level. 



9 http://paperclip. as. arizona.edu/ ~ nsiegler/VLM_binaries 
Maintained by Nick Siegler 
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Only the seeing limited survey of Billeres et al. (2005) found one wide system: Denis 
0551. While it is difficult to estimate the total number of unique VLM objects searched 
we note that Close et al. 2003 & Siegler et al. imaged 80 M6-L0 systems with AO; Reid 
et al. 2006 looked at 52 L dwarfs and Burgasser et al. 2003 looked at 10 T dwarfs with 
HST; Allen et al (2007) has looked at 123 M7-L8 systems. Hence out of ~ 265 unique VLM 
objects (this is a lower limit since not all studies publish their null results), none were in 
wide binaries. Hence we can estimate an upper limit to fvLM mdeilM < 1/265 = 0.37%. The 
lower limit to frequency can be bounded by the fact that one wide system exists, hence 
of the 580 L & T field dwarfs known (Dwarf Archives 0), only one is a wide binary, so 
fvLM md ^ old > 1/580 = 0.17%. Hence we adopt f V LM wMM ~ 0.3 ± 0.1%. This is much 
smaller than our young wide VLM frequency of 6 ± 3%. While a more careful analysis of the 
problem is warranted, it is still clear that the frequency of wide (> 100 AU) VLM binaries 
(Mtat < 0.2M Q ) is ~ 20^ times higher for young (< 10 Myr) VLM objects than for old 
field VLM objects. 

A small observational selection effect is that a much older (~5 Gyr) analog of Oph^llB 
maybe too cool (very late T or "Y" spectral type) and too faint to be easily detected in the 
field today. However, older Oph#16AB, Denis 1618AB, and USCO 1600AB systems in 
the field would be detected by the 2MASS survey as L dwarf primaries with L or early T 
dwarf secondaries (detectable by the aforementioned VLM binary surveys). Correcting for 
this selection effect (throwing out Oph#ll) yields a slightly smaller "field detectable" wide 
VLM binarity of 3/68 = 4.4 ± 2.5% and hence a corrected X evap = 15^^. The lack of old 
wide systems is real, it is not a selection effect of our wide young binary population being 
too low in mass to be detected in the older field. 



4-6.3. Do Wide VLM Binaries Evaporate as they Age? 

One amelioration of the low (~ 0.3 ± 0.1%) wide binary frequency in the field and the 
higher frequency of young (< 10 Myr) wide VLM binaries is that these wide VLM systems 
are dynamically dissipating, or evaporating, over time. After all, a differential "kick" of 
order V esc ~ 1 km/sec would be more than adequate to dissolve all the wide VLM binaries 
in Figure 15. 

Why would such a tidal kick occur? Over the lifetime of a wide binary there will be 
many smaller stochastic encounters (with other stars or molecular clouds; Weinberg et al. 



10 http://spider. ipac.caltech.edu/stafF/davy/ARCHIVE/ 
Maintained by Davy Kirkpatrick, Chris Gelino, & Adam Burgasser 
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1987), which can increase the separation of a wide binary slowly over time. Eventually, these 
encounters may also disrupt the binary. 

However, it is not certain by which mechanism VLM binaries will become unbound. 
Separations of ~ 200 AU are still very small compared to the average 10 5 ~ 6 AU separations 
between stars in the galaxy (whereas wide stellar mass binaries have separations of ~ 10 4 
AU, and so are more easily disrupted). Hence, the odds of any one encounter being close 
enough (< 240 AU) to tidally dissolve a wide (200 AU) VLM binary are very low (Close et 
al. 2003; Burgasser et al. 2003) unless the stellar density is high. 

To investigate the stability of wide binaries we note that Weinberg et al.(1987)'s ana- 
lytic solution of Fokker-Planck (FP) coefficients describing advective diffusion of a binary 
due to stellar encounters is U(a a ) ~ 3.6 x lO^n./O.OSpc" 3 )" 1 * (M tot /M e ) * (MJM Q )- 2 * 
(V re i/20kms~ 1 )(a o / AU)~ 1 Gyr where £*(a D ) is the time required to evaporate a binary of an 
initial semi-major axis of a Q , the number density of stellar perturbers is n* of mass M* 
and relative velocity V re i (adopted from Weinberg (1987); assuming, as they do, that their 
Ink ~ 1). Hence, the maximum projected separation of a bound binary (assuming semi- 
major axis a = 1.26 x sep; Fischer & Marcy 1992) after 10 Gyr in the field is given by: 

Effusive* < 2g 1Q 3 ( J>1*\ 1 ( MuA / 0J V 2 ^ ( J4^\ A (1) 

Ffield ~ \0.05pc~ 3 J \M Q J\M Q J \0.1M e J KJ 

where we have used the measured Galactic disk mass density of O.llM /pc~ 3 and an 
average perturber mass of 0.7M Q , and V re i ~ 20 km/s (Pham et al. 1997; Holmberg & Flynn 
2000). 

In addition to the evaporation of binaries due to diffusion there is also the chance of 
a catastrophic encounter evaporating the binary. While such encounters are less important 
than diffusion, they cannot be completely ignored. From the work of Weinberg et al. we find 
the maximum projected separation (sep) of a binary to stay bound w.r.t. close encounters 
over 10 Gyr in the field is: 

< 52 x 10- (jJgLj) (^) (H) ~ 3300 (J**-) AU (2, 

It is clear that all field VLM binaries with M tot > O.1M will be stable against any type 
of stellar encounter as long as sep < 1800 AU. To better illustrate these regions of stability 
we plot the above two relations on the right-hand side of Figure 17. Note how in Fig. 17 
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all the VLM binaries are not in the "Field Unstable" region (since sep « sep dt /J^ slve * < 
se P C fieid rOPhtC *) ■ Therefore, we can conclude that no known VLM binary will be evaporated 
due to random stellar encounters in the field. However, these limits do help us understand 
the distribution of wide stellar mass binaries like those of Close et al. 1990 (solid triangles 
in Figure 17). Indeed no stellar binaries are observed in the "Field Unstable" region as one 
would expect. 

This still leaves us with the problem of why X evap is ~ 15iiQ. Besides being young, 
all the known wide VLM systems are in the Chamaeleon I (2M1101), or Ophiuchus/Upper 
Sco/Sco-Cen OB (Oph 11, Oph 16, Denis 1618) clusters. Since 2M1101 is near the lower 
core of Chamaeleon I, with ~ 100 members with a core radius of ~ 0.25pc (Luhman 2004) 
we estimate that near 2M1101 ~ 1500/pc 3 . For the Ophiuchus core there are > 200 
members within a radius of ~ 0.3 pc (Wilking et al. 2005) and so n* ~ 1800/pc 3 . At the 
distance of Oph#ll and Oph#16 there are ~ 100 more members but the density drops to 
n* ~ 50 — 220/pc 3 (Wilking et al. 2005). In summary, we adopt a mean ~ 1000/pc 3 
for these "clusters" (in general agreement with other similar sized clusters; Gutermuth et 
al. 2005). These "clusters" (or associations) are much higher densities than the field. It is 
also true that the encounters have much longer interactions timescales (since V re i < 3 km/s), 
hence the clusters where these stars formed will play a role in evaporating them before they 
can join the field. If we assume that these objects are subjected to an additional ~ 10 Myr 
of the mean cluster density then from the above equations we have: 



^//™ < 2g x 1Q 3 ( 1000 V 1 (MuA ( 0/A - 2 (_3_\ „ 44 ( AU 
Cluster ~^ X1U ^ .05pc-3 ) { M Q J \M q ) \2Mm / s ) ^ \0.1M Q ) AU 

(3) 

In the case of a catastrophic encounter we have: 



catastrophic* < „ v , n 3 ( ^ \ ' f M t ot\ ( 0.7 \ 2 ( 3 \ / M tot \ 

(4) 

Hence, after just 10 Myr the cluster has a much more disruptive effect on the wide VLM 
binaries than do encounters in the field over 10 Gyr. From Figure 17 we see that all known 
wide (sep > 100 AU) young (< 10 Myr) VLM binaries are found in the "Cluster Unstable" 
region where sep^l^ e s r tve * < sep < sep d j/J^ stve * . Therefore, there is a good chance that all 
the known wide, young, VLM binaries are in the process of evaporating in their clusters. 
Indeed past cluster member (non-catastrophic) encounters may have already played a role 
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in creating the current (very wide) separations observed for this handful of objects. 

Our simple analysis above is useful but approximate. However, an independent analytic 
solution of this problem by Ivanova et al. (2005) and Brandeker et al. (2006) implies that 
a Mtot < 0.1M Q binary is "soft" if a > 64 AU. Hence, all our wide, young, VLM systems 
can be evaporated by a strong encounter. The timescale for such an encounter by their 
equations is larger than found above (our equation 4), yet their solutions still predict that 
Oph 11 would become unbound in ~ 11 Myr similar to our expectations above. 

Of course none the analysis above implies that all these binaries must evaporate, only 
that it is possible. A more detailed numerical analysis of each system, in each cluster, would 
have to be carried out to ascertain individual outcomes. Indeed some of these may escape 
their cluster before evaporation (as may have been the case for Denis 0551). 

If we assume that ~ 10% of VLM objects form in very low mass groups (n* < 10/pc 3 ), 
then their wide systems will likely survive (whereas none survive for the 90% formed in 
clusters). Then roughly ~ 6 ± 3% of these isolated VLM objects should be in wide systems. 
Hence, one expects ~ 0.6 ± 0.3% of the field to be composed of wide VLM systems. Since 
this is consistent with our observed value of fviM -j , u ~ 0.3 ± 0.1% we can take some 
comfort that the above argument might naturally explain why X evap is observed to be 151^. 

Our scenario does predict that once the system's birth clusters themselves dissolve that 
FyiM wlde ~ FvLM widegcoU - m other words, once the VLM binaries leave their clusters, the 
fraction of wide VLM binaries is "frozen" at the FvLM wid ^ M ~ 0.3 ± 0.1% level. We can 
test this prediction by looking in low density open clusters of intermediate ages to test if 
FvLM wide ~ 0.3 ± 0.1%. The ideal cluster for such a search would be the Pleiades (age ~120 
Myr; d~120 pc). Several large HST & AO surveys have failed to detect a single wide VLM 
binary in this cluster (Martin et al. 2003; Bouy et al. 2006c). A similar lack of wide system 
have been found in the older Hyades (Siegler et al. 2003). So by ~ 100 Myr there is some 
evidence that the wide binary population has mainly been evaporated, as Figure 17 and 
equations 3 & 4 would predict. 

For completeness we note that disruption by molecular clouds does not likely play a 
significant role, since the critical FP impact parameter b'pp 10 ~ 1 x 10 5 AU which is too 
wide to effect ~ 200 AU VLM binaries. However, the much wider (~ 2 x 10 4 AU) stellar 
mass binaries are certainly affected by molecular clouds and their R ~ 2 — 3 pc sub-clumps 
(Weinberg 1987), but they are not relevant in the VLM regime. 
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5. Conclusions 

We have obtained the first spatially resolved R ~ 2000 near-infrared (J & K) spectra, 
mid-IR photometry, and orbital motion estimates for the wide brown dwarf binary Oph #11 
(Oph 1622-2405). We estimate for 11 A and 11B gravities (log(g)> 3.75), ages (5 ± 2 Myr), 
luminosities (log(L/L o )=-2.77±0.10 and -2.96±0.10), and temperatures (T eff = 2375±175 
and 2175±175 K). We find self-consistent DUSTY evolutionary model (Chabrier et al. 2000) 
masses of l7~^M Jup and \^%M Jup , for 11 A and 11B respectively. Our masses are higher 
than the previously reported 13-15 Mj up and 7-8 Mj up masses of Jayawardhana & Ivanov 
(2006b). Hence, we find the system is unlikely a "planetary mass binary", but it has the 
second lowest mass and lowest binding energy of any known binary. 

The #16 binary (Oph 1623-2402) is also an unusually wide (projected separation of 
212 AU), very low-mass (VLM), binary composed of a ~ 100Mj Mp primary (16A) and a 
~ 73M Jup (16B) secondary. 

While young VLM datasets are dominated by small number statistics, they are sugges- 
tive of a moderately common (6 ± 3%) young (< 10 Myr) population of wide (> 100 AU), 
very low-mass (M tot < O.2M ), binary systems like Denis PJ161833.2 -251750.4, 2MASS 
J11011926-7732383, Oph#ll (Oph 1622-2405), and Oph#16 (Oph 1623-2402). Our esti- 
mate of the wide young binary fraction agrees well with the 51^% found in the Upper Sco 
by Bouy et al. (2006). We argue that, over time, the majority of these wide VLM systems 
become unbound leaving just the surviving (most tightly bound) fraction observable in the 
field today. 

We suspect these wide systems likely become unbound since only ~ 0.3 ± 0.1% of field 
VLM objects are wider than 100 AU today. Hence it appears that the field population is 
depleted in wide VLM binaries by a (selection effect corrected) factor of X evap = 15^q- The 
exact mechanism for this evaporation of wide VLM systems is unlikely to be interactions with 
stars in the field or with molecular clouds and their sub-clumps. However, from the Fokker- 
Planck solutions of Weinberg et al. (1987), we find stellar encounters with "birth cluster 
stars" may be efficient in dissolving these wide binaries. Indeed such encounters may help 
boost the observed sep past the sep c ^^ f ophlc * ~ 80 (oTmq) AU limit where catastrophic 
cluster encounters can occur. More detailed numerical studies will be required to better 
constrain the evolution of this evaporating population. 
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Table I. The Oph #11 Brown Dwarf Binary (Separation= 1.943 ± 0.022"; 

PA=176.2±0.5°) 



Parameter 


J 


H 


Ks 


Ls 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


A mag 
11A (mag) a 
11B (mag) a 


0.82 ±0.03 
15.04 ±0.05 
15.86 ± 0.05 


0.69 ±0.03 
14.19 ± 0.07 
14.88 ± 0.05 


0.52 ±0.03 
13.92 ± 0.07 
14.44 ± 0.44 


0.36 ± 0.04 


0.24 ± 0.04 
13.24 ±0.04 
13.48 ± 0.04 


0.11 ±0.04 
13.08 ±0.03 
13.19 ±0.04 


-0.20 ±0.05 
12.96 ±0.05 
12.76 ±0.05 


-0.51 ± 0.10 
12.84 ±0.11 
12.34 ±0.08 



a the J, H, and Ks fluxes are standard MKO magnitudes determined at the Gemini telescope. The Spitzer magnitudes at [3.6], [4.5], [5.8], 
and [8.0] /im are based on the standard Vega IRAC zcropoints of 276.79, 179.5, 116.69, and 63.122 Jy respectively. We determined new IRAC 
photometry with a proper binary PSF fit to the IRAC images. There may be an additional ~15% absolute calibration uncertainty in the 
IRAC photometry (Evans ct al. 2003). 



Table 2. The Oph #16 VLM Binary (Separation=1.696 ± 0.005"; PA=218.29 ± 1.00°) 



Parameter 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] (mjy) 


A mag 
16A (mag) b 
16B (mag) b 


0.691 ± 0.030 
11.19 ± 0.26 
11.94 ±0.26 


0.691 ±0.022 
10.54 ±0.17 
11.23 ±0.17 


0.757 ±0.029 
10.20 ±0.12 
10.88 ±0.12 


0.65 ±0.03 
9.88 ±0.03 
10.53 ± 0.06 


0.67 ±0.03 
9.60 ± 0.03 
10.27 ±0.05 


0.75 ±0.04 
9.10 ±0.03 
9.85 ±0.04 


0.75 ±0.04 
8.12 ±0.16 
9.31 ±0.10 


_ a 

32±5 a 
16±3 a 



a A[24] was not determined, however, since 16B is consistently ~0.75 mag fainter we have estimated that this ratio holds for the[24] ^m 
mjy flux as well. 

b Above magnitudes were dereddened by Ay = 3 ± 1 mag (similar to the Ay = 2 of Allcrs ct al. 2006). The error in the extinction 
dominates the J, H, and Ks photometric errors. With the Ry =3.1 extinction law of Fitzpatrick (1999) we have dereddened our observed 
Oph 16A & B magnitudes by 0.78 mag at J, 0.51 at H and 0.36 at Ks for Ay = 3. Our integrated Oph 16 flux calibration is based on 
integrated flux values from the 2MASS point source catalog and the integrated IRAC values are from Allers et al. (2006). 
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Table 3. The Young, Wide, Low-Mass Ophiuchus Binaries (D=125 ± 25 pc) 



Allers 


2MASS 


SpT 


T <=// 


Lum. 


Gravity a 


Age 


Mass 


Pro. Sep. 


Period b 


# 


System Name 




(K) 


log(L Q ) 


log(g) 


(Myr) 


(Jup) 


(AU) 


(xlO 3 yr) 


11 A 


16222521-2405139 


M9±0.5 


2375 ± 175 


-2.77 ±0.10 


4.25 ±0.50 


5±2 


17.5 ± 2.5 C 


243 ± 55 


20 ± 10 


11B 




M9.5±0.5 


2175 ± 175 


-2.96 ± 0.10 


4.25 ±0.50 




15.5±2.5 C 






16A 


16233609-2402209 


M5±3 d 


3000 ± 300 


-1.18 ± 0.10 




~ 1 


~ 100 c 


212 ± 43 


~ 8 


16B 




M5.5±3 d 


2925 ± 300 


-1.47 ±0.10 






~ 73° 







a Surface gravities estimated based on the best synthestic spectra fits to our J and K spectra. However, a gravity of log(g)~ 3.95 
corresponds to the best fit of Oph#llA and Oph#llB to the spectral standards. 

b Periods estimated based on face-on circular orbits. 

c Masses of 11A and 11B determined from the full range of masses consistent with all of our determined range of gravities, ages, 
luminosities, and temperatures from the DUSTY tracks (Chabrier et al. 2000). No additional error due to unknown systematics in the 
models themselves have been added, hence these mass errors are likely underestimated. 

d Since no spectra have been obtained for 16A and 16B, these spectral types are simply estimated from our dereddened NIR colors from 
Table 2. 

°Masscs of 16A and 16B estimated from the 1 Myr isochrone. Without a detailed spectroscopic observations it is impossible to be more 
accurate about the ages or masses of Oph 16 at this time. 
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OphU Ophl6 




Fig. 1. — Images of both binary systems found from our survey of Ophiuchus VLM objects. 
All images are in the Ks band. The Oph #11 image is seeing limited (FWHM=0.3") with the 
NIRI IR camera (0. 02177" pix" 1 ) at the Gemini North Telescope. Whereas images of Oph 
#16 were obtained at the Keck II telescope with NIRC2 camera (O-OWpix^ 1 ) and the LGS 
AO system. Locking on Oph #16 as its own tip-tilt star gave excellent FWHM = 0.05" 
LGS AO (Wizinowich et al 2006) corrected images. North is up and East left in all images. 
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Fig. 2. — Images of the Oph #11 system at J, H, & Ks (obtained in excellent FWHM ~ 0.3" 
seeing at Ks) with the NIRI IR camera at the Gemini North Telescope. The Ls image is 
from the NIRC IR camera (0.15"pix _1 ) at the Keck I telescope. North is up, east to the left. 
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Fig. 3. — IRAC [3.6] micron (1.2"pix~ l ) image of #11 (upper left; note north is 80° coun- 
terclockwise from up). In the lower left we show the residual (stretch intensified by 5 times 
for clarity) of the system fit to 1 PSF (clearly Oph #11 is not well fit by a single PSF). The 
second column from the left shows our #11 binary model (Sep = 1.94", PA = 176°) with 
the best fit fluxes (in this case with a A [3.6] of 0.24 mag). The next (3 rd ) column from the 
left shows a single star and its residual when fit with our IRAC PSF (giving a similarly good 
residual as our binary model). The upper image in the right-hand column shows a simulated 
system similar to the #11 AB binary and then below it the residuals left after fitting a single 
PSF fit (note the similarity of this simulation to the comparable fit to the actual data shown 
in the left-hand column). The robustness of the binary fit gives confidence in our IRAC 
photometry to the ~0. 03-0. 05 mag level. Similar fits were also obtained at [4.5], [5.8], and 
[8.0] p,m for #11 and #16; see Tables 1 & 2 for the resulting mid-IR photometry. 
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Fig. 4. — The J-H, H-Ks color-color diagram for field (old) low-mass M, L, & T dwarfs. Note 
how the positions of 11A and 11B are similar to that of AB Pic b another very low- mass, 
young, object (Chauvin et al. 2005b). The offset in the positions of 11A and 11B from the 
dwarf locus is likely due to the low surface gravity of these objects compared to field dwarfs 
(Song et al. 2006). Figure adopted from Song et al. (2006). 
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Wovelength [yU.m] 

Fig. 5. — The K-band NIRSPEC spectrum of Oph 11A (red line). Note the good agreement 
between the young LO dwarf 2M0141 (from Kirkpatrick et al. 2006) and 11A. In particular, 
the good fit to the temperature sensitive CO lines suggests 11A is just slightly warmer than 
a young LO. The poorer fit of the gravity sensitive Nal and Cal to the older (30 Myr) star 
GSC8047 (from Chauvin et al. 2005a) suggests that 11A must be younger than 30 Myr (we 
have lowered the resolution of the Oph 11A spectra to match that of GSC8047). 
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Fig. 6. — Similar to Figure 5, except for 11B. Note the even better fit to 2M0141 for 11B. 
Hence we suggest that 11B is similar to a young M9.5 ± 0.5, consistent with the optical 
spectral type of M9.5-L0 found by Jayawardhana & Ivanov (2006b). 
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Fig. 7. — The J-band NIRSPEC spectrum of Oph 11A compared to a range of young brown 
dwarfs. Here we see the many gravity sensitive features in the J-band (McGovern et al. 
2004) can be used to estimate the gravity (and hence age) of 11A. To minimize uncertainly 
we compare to other spectra obtained with the same instrumental set-up as we used. These 
J-band comparison spectra from the work of McGovern et al. (2004) and Kirkpatrick et al. 
(2006) show that llA's KI doublets and J-band pseudo-continuum best fit a brown dwarf of 
age 5 ± 2 Myr such as a Ori 51 (McGovern et al. 2004). Moreover, the strength of llA's KI 
doublets (and poor fit to the pseudo-continuum) compared to the very young brown dwarf 
KPNO Tau-4 strongly suggests that the system age for Oph 11 is > 2 Myr. 
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Wavelength [fim] 

Fig. 8. — Similar to Figure 7 but with 11B. We see the surface gravity of the 5 ± 2 Myr 
a Ori 51 is the best fit. Moreover, the poor fit to KPNO Tau-4 strongly suggests that the 
system age for Oph 11 is > 2 Myr. 
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Fig. 9. — Comparison of our observed K-band spectra (red) to synthetic spectra (black) 
computed using up-to-date PHOENIX DUSTY atmosphere models (Hauschildt et al. in prep). 
See Section 3.3.1 for a discussion. 
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Fig. 10. — Same as Figure 9 except for the J-band instead of the K-band. See Section 3.3.2 
for a discussion. 
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Fig. 11.- Synthetic 2MASS and IRAC fluxes for log(g)=3.95 and T eff = 2375K and 
T e ff = 2175K. These SEDs were computed using up-to-date PHOENIX DUSTY atmosphere 
models (Hauschildt et al. in prep), assuming 125 pc as the distance to Oph #11. These 
SEDs demonstrate the significant IR excess for 11 A (red points) and 11B (blue points) at 
A > 3/im. Component 11B appears to have a much stronger excess than 11A. The integrated 
24/zm datapoint for the system is plotted since we could not resolve the system at 24/xm. 
Overall, the strong warm dust excess and Ha emission suggests ages < 8 Myr, consistent 
with our estimated 5 ± 2 Myr age. 



-40- 




3000 ! ^ 2500 ~ 2000 

Teff (K) 



Fig. 12. — An HR diagram where the dashed "vertical" lines are iso-mass contours for 
the DUSTY models (from left to right: 94, 73, 52, 42, 21, 16, 13, 10, & 7 M Jup ), while 
the four slightly more "horizontal", lines are the DUSTY isochrones (top to bottom: 1 
Myr (dash-dot) and 5 Myr (solid line), lOMyr (dash-dot) and 100 Myr (dash-dot)). The 
2MASS 0535 data points (open circles) mark the only known low-mass, very young, binary 
with a well determined dynamical mass. The thick short diagonal lines representing the 
displacement from the measured luminosity and temperature to the values actually predicted 
by the DUSTY models. One can see that for the secondary of the 2MASS 0535 eclipsing 
binary system (in Orion; Stassun et al. 2006) the models have some error, but in general the 
binary is close to the relevant 1 Myr isochrone. Hence we can have some confidence in these 
models to estimate masses. From the intersection of the 11 A (blue area) and 11B (red area) 
areas w.r.t. the estimated age of 3 — 7 Myr (from our NIR surface gravity measurements) we 
estimate masses of 13 — 21Mj up and 10 — 20Mj up for 11A and 11B. Therefore, the system 
appears to be a very low-mass brown dwarf binary. 
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Fig. 13. — DUSTY tracks of Chabrier et al. (2002) in terms of surface gravity (log(g) in 
cgs units) against temperature (T e ff). Hence, we can estimate masses independent of the 
possibly variable NIR luminosity of Oph 11. From our synthetic spectral fits (with the same 
models) we find log(g)> 3.75. Moreover, the existence of a strong IR excess and Ha emission 
(Jayawardhana & Ivanov (2006b)) suggests that the age of Oph 11 is < 8Myr. Therefore, 
we "shade-in" the remaining areas of model space that are self-consistent for 11A and 11B. 
For 11A (red "area") we find 17~^Mj up is the allowed mass range. For 11B (blue "area") we 
find 14^lMj up is the allowed mass range. 
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Fig. 14.— Upper left a 1993.33 epoch red POSS II image of Oph 11A and Oph 11B (11B 
is brighter than 11A in this filter due to its stronger Ha emission). Note how nearby field 
PSF stars have circular PSFs very different from the Oph 11 system. This image shows the 
system with a separation of 2.05 ± 0.20" and a PA = 174 ± 6° which is consistent with the 
current (2006.55) orientation and separation of the binary. Hence, we measure only ~ 3 ± 5 
km/s of "motion" in the plane of the sky over the last 13.22 years, consistent with the system 
being a bound, common proper motion, pair. This image has been unsharp masked, and 
magnified 5x5. The bottom panels show the IR POSS II image of the Oph 16 system in 
1982.66. We find in 1982.66 the system had a separation of 1.87 ±0.2" and PA= 216.3 ±10°, 
this position is consistent with its current position 23.89 years later. This implies the Oph 
16 system is likely also a bound common proper motion system. North is up and east is left 
in all images. 
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Fig. 15. — The Oph 11 system is compared to other known VLM and brown dwarf binaries 
(old VLM systems are open stars, open circles are young (< 10 Myr) VLM systems, and 
stellar mass binaries are solid symbols). Adopted from Close et al. (2003) and Burgasser et 
al. (2006). 
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Fig. 16. — Illustrating the uniquely low binding energy of the Oph #11 system compared to 
other known VLM systems. Note how our newly discovered Oph #16 binary is also weakly 
bound. Adopted from Close et al. (2003) and Burgasser et al. (2006). 
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Fig. 17. — Same as Fig. 15, but with the instability zones highlighted. The zones are 
determined by equations 1-4 which predict the maximum bound separations plotted as 
dashed "diagonal" lines from left to right as sepj^™*, sepf^ r oph ' lc * sepf/J™* , and 
se P C fieid r ° phlC * ■ Note how all young (< 10 Myr; open circles), wide, VLM systems are in 
the "Cluster Unstable" region. Therefore, we expect that many of these "cluster" (n* ~ 
1000/pc 3 ) binaries could be evaporated before joining the field (open stars > 100 AU). 
Hence, wide VLM binaries in the field should be rare -as observed. We suggest that only 
wide VLM systems in very low (n* < 100/pc 3 ) density groups will survive and join the field. 



